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Abstract 

By considering the local vector current between nucleon states and im- 
posing charge conservation, we determine its renormalisation constant and 
quark mass improvement coefficient for Symanzik 0(a) improved Wilson 
fermions. The computation is first performed for quenched fermions (and 
for completeness also with unimproved fermions) and compared against 
known results. The two-flavour unquenched case is then considered. 
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1 Introduction 



A naive discretisation of fermions onto a hypercubical lattice gives for the action 
errors of 0(a 2 ). However returning to the continuum then leads to the famous 
fermion doubling problem when we find 15 extra copies of our original fermion. 
To cure this Wilson [l] added the 'Wilson term' to the action so that the copies 
decouple in the continuum limit: but then discretisation errors are 0(a). As the 
gluon part of the action (sum of plaquettes in this article) has already 0(a 2 ) 
errors, it is desirable to also achieve this for the fermion action. The Symanzik 
programme 1 allows a systematic reduction of errors 2 to 0(a 2 ). An additional 
operator of dimension 5 (the 'clover term') is added to the Lagrangian with a 
coefficient c sw suitably adjusted so that on-shell quantities such as masses now 
have 0(a 2 ) errors. For matrix elements it is also necessary to add further higher 
dimensional operators to the original operator to achieve 0(a) improvement. In 
this letter we shall be concerned with the determination of the improvement 
coefficients of the associated improvement operators for the local vector current: 
yw = (pf^q. For this operator just two additional operators am q V^ and iadyTjfJ 
are required giving the 0(a) improved vector current yw /MP , and renormalised 
current Vy' R , as 

V<«> fl = Z V V^ IMP , V^"" = (1 + am q b v )(V^ + iacyd v T$) , (1) 

with Tj$ = qa^q, o> = «[7/*,7«/]/2, and d^(p(x) = [<f>(x + fi) - <f){x - p)]/(2a). 
The iad v T^} operator only plays a role in non-forward matrix elements and will 
not be considered further here. Thus to 0(a) improve the local vector current 
we need to determine the improvement coefficient bv(go) (where go is the bare 
coupling constant). Also as this current is not conserved on the lattice then, as 
discussed in the next section, it is renormalised with renormalisation constant 
Zy(go)- Perturbatively we have El [HI| to one loop (independently of the 
presence of fermions), 

Z v (g ) = 1 - (0.174078 - 0.040069c sw - 0.004586c;; J # 2 + . . . , (2) 

where for unimproved fermions c sw = and for 0(a) improved fermions c sw = 
1 + O(<7o), together with 

b v (g ) = l + 0.15324g 2 + ... . (3) 

However in regions where numerical simulations are performed, (3 = Q/g^ ~ 
6.0 — 6.4 for quenched fermions (nj = 0) and 5.2 — 5.3 for unquenched fermions 
(rif = 2), the above formulae may not be applicable. 

1 For an introduction see, for example, 0E1E]- 

2 Ginsparg- Wilson fermions, an alternative formulation, automatically have 0(a 2 ) errors 
|H] and have better chiral properties but are very CPU time consuming. 
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In this article we shall determine Zy and by non-perturbatively by considering 
(nucleon) matrix elements of the time component of the local vector current and 
imposing charge conservation (the details will be given in the following section). 
We shall first consider quenched fermions and compare our results with known 
results in the literature. (While most of the results will be for 0(a) improved 
fermions, for comparison, we shall also briefly consider the unimproved case, 
c sw = 0.) This is then followed by the unquenched case. Preliminary results have 
appeared in [TTj . 



2 The Conserved and Local Vector Currents 

There is an exact global symmetry of the lattice action q — > e~ %OLq q, q — *> e lotq q. 
This global symmetry is flavour conservation (if you just rotate the quarks of one 
flavour) or baryon number conservation (if you rotate all quark flavours equally). 
Separate quark transformations are possible because in pure QCD there are no 
flavour changing currents. Upon using the Noether theorem this symmetry gives 
an exactly conserved vector current or CVC of 

J^\ x + IP) = 3 [&(7/i - l )U^(x)q x+(l - g x+A (7/, + 1)^(2;)^] . (4) 

(Being conserved this current requires no renormalisation constant and is 0(a) 
improved 3 .) By this we mean that the Ward Identity, WI, is 

+ m,» - (^) iv + {^) qv , (5) 

where fl is an arbitrary functional of the U, q and q fields. is the backward 
derivative, A fl (p(x) = [4>(x) — <p(x — fi]/a. (Although in eq. © we integrate 
out the fermion fields, (. . .) q , and take the average over the gauge fields, (■ ■ -)u, 
the equation is already true configuration-by-configuration.) We immediately see 
that if the region over which Q is defined does not contain x the RHS of eq. (0) 
vanishes so that A^J^(x + \fi) =0. If Vt contains x then the RHS effectively 
'counts' the number of quarks and anti-quarks in Q. 

Here, in this study, we shall take Q — > B(t)B(0) where B is the standard 
(stationary, or p = 0) nucleon operator, containing two u quarks and one d quark 
summed over the spatial planes 4 . (This was previously used as part of a project to 
determine moments of structure functions; other choices are of course possible.) 



3 This is of course only true when considering forward matrix elements. For non- 
forward matrix elements the CVC requires the additional operator j/ t ^ — > = + 



d u TjfJ(x)+d l/ TlfJ(x + fl) 



4 We actually take T u ^ po1 B (t)B a (O) where T un P° l = |(1 + 74) projects out the unpolarised 
component of the nucleon field, B a (t) = ^-e yfc u^(af,t)[u : '(£,t) T C75d' c (af,i)]. 
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B B B B 




Figure 1: The quark-line connected diagram, left hand picture, and quark-line dis- 
connected diagram, right hand picture. The cross denotes the operator insertion 0(t). 

Inserting this Q in eq. (0) and also summing over spatial points of x = (x, r) 
gives 

R[4 q) (r)] - R[4 q) (r - 1)] = - X W (<5 T , t - 5 Ti0 ) , (6) 
where R is denned as the ratio of three-point to two-point correlation functions, 

and y (u) = 2, = 1. With 



or r(0)/y\i _ J C S ?) 0<r<t 



for an N$ x N T lattice with cf' (i — 1, 2) being constant eq. © can be solved to 
give the result 

c[ q) - c {q) = x {q) ■ (9) 

So R should be a constant with jump, or discontinuity, given by x (q) - Note 
that this should be true to 'machine accuracy'. (Indeed in this special case R 
may also be taken to be the ratio of three- to two-point correlators for a single 
configuration.) This result may also be shown using transfer matrix methods, as 
indicated in the appendix, where it is also demonstrated that c{ ■ 

The lattice computation of the three- and two-point functions for R[0(t)} 
follows the standard way, see eg [T2"] . (The source B(0) and sink B{t) have 
been additionally improved by non-relativistic projection and Jacobi smearing to 
increase the overlap with the ground state nucleon, as described for example in 
[13 . This does not affect the arguments given above.) We must, in principle, 
compute a quark-line connected contribution (to the operator) and a quark-line 
disconnected term as shown in Fig. ^ This latter term is numerically extremely 
difficult to compute, due to ultra-violet fluctuations. However for the CVC this 
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term is in fact constant. This may be easily seen by substituting Q — > (B(t)B(0)) q 
into the WI, eq. ©• The RHS is then zero and so using the same argument as 
before the appropriate ratio is constant for all r. There is thus no contribution 
to the discontinuity. Indeed on a finite lattice, we would expect this constant to 
be exponentially small (with exponent proportional to Nt)- Physically there is 
no quark-line disconnected term because creating a quark- ant iquark pair cannot 
change the charge. Nevertheless with an eye on the computation of the local 
current it is useful to consider the difference between the u and d operators, ie 
the non-singlet operator, 

= q(u) _ (d) ^ ^Q) 

in which the quark-line-disconnected terms cancel and so we find 

D[7 M)/ NT { Cl'^ 0<T<t ( s 

{ c 2 t < r < Nt — 1 

where the discontinuity between the two constants, Ai^fJ^""^], in eq. (fTTj) is 
given by 

AR[4 u ' d) ] = c { r d) - c^ d) = x [u) ~ X (d) = 1 • (12) 

An example for the ratio R[J < £ l ~ d ' 1 ] is shown in Fig. El For jj" - ^ a very good signal 
is observed. For the jump, as expected, we find 1 to within machine precision. 

The local vector current (LVC) does not obey the WI given in eq. but as 
Jja) = + 0(a) there is an additional term formally of 0(a) on the RHS of 
this equation. However perturbatively expanding eq. (j3J) gives loop graphs with 
ultra-violet divergences ~ 1/a and so this additional term gives a contribution. 
Thus to obey the WI to 0(a 2 ), must be renormalised with renormalisation 
constant Zy and the quark-line disconnected term in Fig. ^may give an 0(a) 
contribution. Again to avoid computing this term we consider the non-singlet 
operator 5 . The necessity for a renormalisation constant is also reflected in the 
fact that A-RfV^" ^] will not be equal to one. This is illustrated in the RH picture 
of Fig. |21 So we can define the renormalisation and improvement constants (Zy 
and by respectively) by demanding that the local current has the same behaviour 
as the conserved current, ie 6 

Z v (l + am q b v ) = (ARIV} 11 ^})' 1 . (13) 



5 Tests from computing R\V^ (r)] for q — u, d separately showed indirectly that the contri- 
bution from the quark-line disconnected term must numerically be very small. So the difference 
between the singlet and non-singlet renormalisation constant must also be very small. 

6 It is more precise to define 

Z v = lim (ARiv}'^})- 1 , 



9 , , „ „rrr(u-dh 



lim — -ln(Ai?[V 4 v 

m g ->o o(am q ) 
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Figure 2: R[j\ (r)] and R[V^ (r)] plotted against the operator position r for 
the quenched (n/ = 0) data set /3 = 6.0, k = 0.1342 on an iVj x iV T = 16 3 x 32 lattice 

with i = 13. Typical fit intervals for V4 d ^ are given by the pairs of vertical dashed 
lines. The relationship between am q and k is given in eq. Q14j). 

Thus upon plotting the data for (AR[v} u against am,, we see that the 

intercept gives Zy while the gradient gives Zyby. Note that Zy has potential 
0(a 2 ) differences to other definitions of the renormalisation constant while by 
has possible 0(a) differences. 

While this procedure is correct in the quenched case, for the unquenched 
case for complete 0(a) cancellation, g in Zy should be replaced by g where we 
have #o = 9o(^ + b g am q ), d4. Thus the LHS of eq. (|T3j) should now become 
Zy(l + am q (by + ^bgj^j^-)). So while the intercept still gives Zy, the gradient 
and hence by is modified. b g is only known to one loop perturbation theory, 
b g = 0.01200n/go + 0(g$), [Hj. Estimating dZy/dgo for unquenched fermions 
from the Pade fit, eq. (fTTj) and Table |U to be ~ —1.3 gives roughly for the 
additonal term a decrease of ~ 0.02. As numerically by will turn out to be ~ 2, 
then this could give a 1 — 2% correction. At present, due to uncertainties in 
estimating this extra term, we shall ignore this small correction factor. 
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^ SW 




N% x N T 


^ configs. 


(AR\V} u ~ d) })- 1 


6.0 


1.769 


0.1320 


16 3 x 32 


0(450) 


8839(41 


6.0 


1.769 


0.1324 


16 3 x 32 


0(560) 


0.8717(3) 


6.0 


1.769 


0.1333 


16 3 x 32 


0(560) 


0.8418(4) 


6.0 


1.769 


0.1338 


16 3 x 32 


0(520) 


0.8240(6) 


6.0 


1.769 


0.1342 


16 3 x 32 


0(740) 


0.8124(8) 


6.2 


1.614 


0.1333 


24 3 x 48 


0(3001 


8691(21 


6.2 


1.614 


0.1339 


24 3 x 48 


0(300) 


0.8503(2) 


6.2 


1.614 


0.1344 


24 3 x 48 


0(300) 


0.8342(2) 


6.2 


1.614 


0.1349 


24 3 x 48 


0(470) 


0.8186(2) 


6.4 


1.526 


0.1338 


32 3 x 48 


0(220) 


0.8624(1) 


6.4 


1.526 


0.1342 


32 3 x 48 


0(120) 


0.8504(2) 


6.4 


1.526 


0.1346 


32 3 x 48 


0(220) 


0.8376(1) 


6.4 


1.526 


0.1350 


32 3 x 48 


0(320) 


0.8253(1) 


6.4 


1.526 


0.1353 


32 3 x 64 


0(260) 


0.8163(3) 



Table 1: Parameter values used in the quenched simulations with improved fermions, 
together with (AR[v£ u ^] ) _1 . The fit intervals chosen were [nj , n u ] and [Nt — n\ , Nt — 
n u ] where [ni,n u ] = [9, 13], [6, 11], [7, 16] for /3 = 6.0, 6.2 and 6.4 respectively. 

3 Simulation parameters and raw results for 0(a) 
improved fermions 

We have made runs at (3 = 6.0, 6.2 and 6.4 for quenched fermions and (3 = 5.20, 
5.25 and 5.29 for unquenched fermions. The run parameters and results for 
(A J R(V A 4 (u_d) ])- 1 are given in Tables □ and El 
The bare quark mass is defined as 

am q (K,g Q ) = - ( - j , (14) 

2 \k K c (g ) J 

where k is the hopping parameter of the simulation. It is necessary to find n c (go), 
ie the critical point where the (bare) quark mass vanishes. From PCAC we know 
that the quark mass is oc rn^ s and we can use this to determine where the quark 
mass vanishes. However a more precise/stable determination was often possible 
if the PCAC quark mass m^ WI was used, so we fitted the dimensionless quantity, 

X = F{r m I q MP ) 

= B 1 r m I q MP + B 2 (r m I g MP ) 2 + ---, (15) 

where X = r^m^ 1 and m IMP = m q (l+b m am q ) is the O(a) improved quark mass. 
However determining am q WI requires a knowledge of the ca improvement coef- 
flcent. While this is known for quenched fermions jT^j for unquenched fermions 
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9 ni 71 


PI 1 ^A9 


1 fi 3 y ^9 

ID A OZi 


<JUUU 


\<l v u or 




5.20 


2.0171 


0.1350 


16 3 x 32 


8000 


UKOCD 


7731(10) 


5.20 


2.0171 


0.1355 


16 3 x 32 


8000 


UKQCD 


0.7537(16) 


5.25 


1.9603 


0.1346 


16 3 x 32 


2000 


QCDSF 


0.8019(08) 


5.25 


1.9603 


0.1352 


16 3 x 32 


8000 


UKQCD 


0.7781(08) 


5.25 


1.9603 


0.13575 


24 3 x 48 


2000 


QCDSF 


0.7560(06) 


5.29 


1.9192 


0.1340 


16 3 x 32 


4000 


UKQCD 


0.8328(07) 


5.29 


1.9192 


0.1350 


16 3 x 32 


5000 


QCDSF 


0.7948(03) 


5.29 


1.9192 


0.1355 


24 3 x 48 


2000 


QCDSF 


0.7747(04) 



Table 2: Data sets used in the unquenched, n/ = 2, simulations together with 
{ARlvt^})- 1 . The fit intervals chosen were [9, 13] and [N T -9,N T - 13]. 



p 




Z v 


Zyby 


by 


6.0 
6.2 
6.4 


0.135201(9) 
0.135803(3) 
0.135744(1) 


0.7799(7) 
0.7907(3) 
0.8027(2) 


1.168(10) 
1.135(06) 
1.116(04) 


1.497(13) 
1.436(08) 
1.391(05) 


5.20 
5.25 
5.29 


0.136072(10) 
0.136287(09) 
0.136368(09) 


0.7304(18) 
0.7349(09) 
0.7420(07) 


1.472(44) 
1.460(31) 
1.411(19) 


2.015(61) 
1.987(43) 
1.902(25) 



Table 3: K c (go), the intercept Zy, gradient Zyby and by for the quenched data sets 
(/? = 6.0, 6.2 and 6.4) and for the unquenched data sets {(3 = 5.20, 5.25 and 5.29). 

this more sensitively affects the value of am^ WI . Thus in this case we have set 
X = (r m ps ) 2 . r in eq. (fTK|) is the 'force' scale, (THj. While r could be omit- 
ted for quenched fermions, for unquenched fermions tq becomes quark mass as 
well as coupling constant dependent and we use the results given in jTTj, supple- 
mented by [THJ. For quenched fermions, we fitted both B\ and B 2 coefficients, 
which practically meant that no b m coefficient was needed, being absorbed into 
the B 2 coefficient, while for unquenched fermions, as we had only three quark 
mass values we set B 2 = and used a tadpole improved value of b m , following 
the prescription given in [TU] (and using the plaquette values given in [H]). This 
gave the results for k c in Table El For quenched fermions these numbers are in 
good agreement with those given in |15j . 

As discussed in section E] we now measure the intercept and gradient of 
(ARivt^})- 1 against am q . The results are shown in Fig. El for both the 
quenched and unquenched data. Very good linearity is seen in the results, which 
enables a precise estimation of the intercept and slope. The results are also given 
in Table El While these extrapolations have no problems the determination of k c 
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Figure 3: Linear extrapolations for (AR[\Q for the quenched data at (3 = 6.0, 

6.2 and 6.4 (upper set of lines, from (3 = 6.4, upper line, to 6.0, lower line) and for the 
unquenched data at j3 = 5.20, 5.25 and 5.29 (lower set of lines, with the highest being 
for (3 = 5.29). 



may be more problematic. To give an estimate of possible further errors, if k c is 
varied by Sk c , then we have 

5Z V ~ byZ v 5{am q ) , 5b v ~ by5(am q ) , with 5(am g ) ~ tt^^^c • (16) 

For example, taking 5n c ~ 0.2 x 10 -4 gives 5(am q ) ~ 0.0005 and 5Z V ~ 0.0006, 
5b v ~ 0.001. These are of the same order or less than the given statistical error. 



4 Zy and by for 0(a) improved fermions 

We are now in a position to give our final results for Zy and by. From eq. (JEJ), 
dividing the gradient in Fig. El by the intercept gives in Table 01 in the third and 
fifth columns Zy and by respectively. 

As for both quenched and unquenched fermions we have three (3 values, we 
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Pi 


P2 


Qi 




Quenched (n/ = 


0) 


Zy 


-0.634 


0.0196 


-0.504 


by 


-0.627 


-0.0444 


-0.781 


Unquenched (rif — 


= 2) 


Zy 


-0.796 


0.0652 


-0.667 


by 


-0.614 


-0.0432 


-0.767 



Table 4: The Pade fit results as defined in eq. (|17|). 



can attempt to make a Pade-type fit of the form 

pboH 1+ ff +r g ° . (it) 

1 + Qi9o 

for Z v (g ), by (go) constrained to reproduce the weak coupling results, eqs. (J2J), 
(J2J). There are thus 2 free parameters. In Table 0] we give the results of the fits. 
Note however that even though we have matched to the weak coupling results, as 
the range of j3 where the numerical results lie is rather small and far away from 
this region, intermediate regions may not be represented so well. We estimate 
total errors on these Pade results to be about \% for Zy for both the quenched and 
unquenched cases and for by, 1% and 2% for quenched and unquenched fermions 
respectively. (Note also that for the unquenched results there is a further error 
in by of 1 — 2% due to an uncertainty in b g as discussed previously in section [21 ) 
Alternative non-perturbative determinations for quenched 0(a) improved fer- 
mions have been given by the ALPHA Collaboration, using the Schrodinger func- 
tional method, the LANL Collaboration, [20 , 21] using other Ward identities 
and the SPQcdR Collaboration [22| using the RI-MOM scheme, [23]. All these 
other methods have reasonable agreement with the results given here, as can 
be seen in the following pictures. (One should remember that Zy definitions in 
particular RI-MOM can vary by at least 0(a 2 ) and by definitions can vary by 
0(a).) In Figs. EJ El we plot our results for Zy and by respectively for both 
quenched and unquenched fermions. As, in particular for unquenched fermions 
the available /3-range is rather narrow and far from the continuum limit the Pade 
extrapolation should be treated with some caution. However it is encouraging to 
note that the ordering of the points is correct. Thus the Pade results should be 
regarded mainly as an interpolating formula between the various (3 values. Note 
however in the quenched case that the Pade results track quite well the ALPHA 
results. 

Some of our earlier results for Zy for quenched fermions can be found in |24j . 
Note also an example of the large increase in the statistical fluctuations shown 
there when matrix elements with moving nucleons (p ^ 0) are considered. 
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Figure 4: Zy (LVC, filled squares) determined in this work for quenched and un- 
quenched 0(a) improved fermions. For quenched fermions a comparison is made with 
ALPHA, QH, LANL, [H] and SPQcdR, [H]. Pade fits are also given for our (full 
lines) and the ALPHA (dot-dashed line) results. The dashed line is the first order 
perturbative result, eq. (J2J). 

5 Zy for unimproved fermions 

Although most of our results are for 0(a) improved fermions we also have results 
at P = 6.0 for quenched unimproved fermions (c sw = 0) which serve as a useful 
check on the improved results. The method is identical to that described previ- 
ously, so we shall just give the results here. In Table El we give the parameters of 
the runs and raw results. Note that the range of quark masses is greater than in 
the 0(a) improved case and that we have runs on different volumes at the same 
quark mass. It is apparent from the table that finite volume effects are rather 
small. 

Plotting (ARlvt^})- 1 against am q gives the result Zy = 0.6436(3) at 
(3 = 6.0 (for completeness the gradient is 0.9088(27)), using k c = 0.157129(10). 
We see that the numerical value is somewhat lower than the equivalent 0(a) 
improved value. This result is consistent, but a little larger than the recent de- 
termination given in j2H] (although there different /5-values are used) . Again this 
is not unexpected, as different determinations of Zy can now differ by 0(a) terms 
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P 

Figure 5: by (LVC, filled squares) determined in this work for both quenched and 
unquenched 0(a) improved fermions. Also shown are our Pade fits (full lines) and that 
from ALPHA, |19j (dot-dashed line), and the LANL, [2^ results for quenched fermions. 
The dashed line is the first order perturbative result, eq. (0). 

(see also [24J). 

6 Conclusions 

Our method is not in disagreement with results of other approaches for 0(a) 
improved quenched fermions. For example even for by there is only a few percent 
scatter in comparison with other methods. This is what one would expect with 
0(a) discrepancy between the various definitions. 

Zy and by are both further away from 1 in the unquenched case than in the 
quenched case at comparable lattice spacings (roughly a n/=2 (5.25) ~ a n/=0 (6.0)). 
This is partially an effect from the use of (numerically) smaller j3 factors and also 
partially because of the presence of fermions which induces a shift in the effective 
(3. Finally we note that we cannot use one loop perturbation expansion results in 
the present day quenched/unquenched (3 regions, although this can be improved 
by using TI-RGI-BPT, eg [H], or BLM/TI methods, £7]. 
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(IAU[V A \) 


n 

U .\J 


1 51 ^ 

U.lul J 


Ifi 3 x 39 




7591 (9) 


6.0 


0.1530 


16 3 x 32 


0(1130) 


0.7220(2) 


6.0 


0.1550 


16 3 x 32 


0(1360) 


0.6845(3) 


6.0 


0.1550 


24 3 x 32 


0(220) 


0.6837(3) 


6.0 


0.1558 


24 3 x 32 


0(220) 


0.6689(4) 


6.0 


0.1563 


24 3 x 32 


0(220) 


0.6595(6) 


6.0 


0.1563 


32 3 x 48 


0(250) 


0.6599(4) 


6.0 


0.1566 


32 3 x 48 


0(410) 


0.6545(5) 



Table 5: Parameter values used in the quenched simulations with unimproved 
fermions, together with (AR[vj u 
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Appendix 

The result in eq. may alternatively be shown using transfer matrix methods. 
As J4 is conserved then (with normalisation = 5^) 

(^i 9) |j) = xl%, (18) 

where xf is the 'charge' (associated with the current J$) of the state Thus 
inserting complete sets of states into the three-point function gives 

(B(t)4 q) (r)B(0)) = Y,(mi)x?\i\me~ aEit e- aEj{NT - t) 0<r<t 

ij 

= (0|B|A0xi?(Af|5|0)e- amiV ' + ••• , (19) 
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(where the second equation assumes that the excited states have died away and 
takes the charge on the nucleon, Xn\ to De Xn = (^l^i^l^O) and similarly 

(4 q \r)B(t)B(0)) = Y,xf(m^(i\B\j)e- aEit e^ N ^ t < r < N T 



1.1 



= xS(^1^|0)(0| J B|A^*)e- am ' v * (iVT -' ) + ... , (20) 

N* being the parity partner of the nucleon and where Xn* = (N*\J^\N*) = 
—X$- All of these expressions are independent of r. The transfer matrix ap- 
proach thus again predicts perfect plateaus for J^ g \ Furthermore, the disconti- 
nuity is given by the difference between these two expressions. As (i\B\j) is only 
non-zero if the difference in the charges of \i) and \j) is equal to that of the charge 
of a nucleon, then we can replace xf^ ~ X^ by X*n to give 

(B(t)J^(r)B(0) )\ 0<r<t - (J^(r)B(t)B(0) )\ t<T<Nr 
= E(xl ?) - X?)(mi)(i\me- aEit e- aE ^ NT ^ 

ij 

= x®(B(t)B(0)). (21) 

From eqs. (fT^j) . (|2T?j) we see that for t < N T /2 (the case considered here) the 
exponential factor is larger for the term < r < t which means that 4 9) is small 
in comparison with c[ . 
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